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Abstract We continue the investigation of how to use the divergence-free con¬ 
dition to resolve the az imuthal ambiguity present in vector magnetogram data . 
In p revious articles, bv ICrouch. Barnes, and Lekal {Solar Phys. 260, 271, 2009l l 
and [C roiichl {Solar Phvs. 282. 107. l2013l l. all methods used an expression for 
the divergence of the magnetic field that involves differentiation of quantities 
that depend on the choice of azimuthal angle. As a result, all heights used to 
approximate line-of-sight derivatives should generally be disambiguated simul¬ 
taneously. In this article, we investigate a set of methods that use an expression 
for the divergence that involves differentiation of quantities that do not depend 
on the choice of azimuthal angle. This results in an expression for the diver¬ 
gence that can be used to disambiguate each height independently. We test two 
methods using synthetic and find that the two-step, hybrid method, adapted to 
disambiguate each height independently, generally produces reasonable results. 
Moreover, the time required to compute solutions is substantially decreased in 
comparison to the corresponding method that disambiguates all relevant heights 
simultaneously. 


Keywords: Sun: magnetic field 


1. Introduction 


Reliable measurements of the physical conditions in the solar atmosphere, such 
as the magnetic field, are required to constrain models for a wide range of solar 
phenomena. To fully determine the magnetic field vector, however, involves the 
resolution of the azimuthal ambiguity, which arises when the component of the 
field in the direction perpendicular to the line-of-sight is inf erred from the linear 
polarisation of magnetically sensitive spectral lines {e.g., Harvey, 1969h . This 
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approach cannot distinguish between the two possible choices for the perpen¬ 
dicular component of the held, which have the same magnitude but differ in 
direction by 180°. 

Several methods are currently in use for resolving the azimuthal ambiguity 
in single-height vector magnetog ram data; a selection of t hese were r e viewe d 
and tested on synthetic data by Metcalf el al\ ( 200(jl and Leka et ai\ ( 2009h . 
Single-height vector magnetogram data include information about the variation 
of the magnetic held over the solar surface, but not in the direction perpendicular 
to it. Therefore, disambiguation algorithms that need this information (say, to 
compute the divergence of the held) must make an assumption. For example, 
some versions of the minimum energy method use a potential-held extrapolation 
to approximate the variation of the magnetic held in the direction perp endicular 
to the solar surface, which is needed to compute the divergence (e.q.. Metcali 
1994 : Metcalf et ai . 20061 : Leka. Barnes, and Crouch . 20091 : Leka et all l2009fl . 


This approximation may not be an accurate representation of the solar magnetic 
held in general. 

Another class of methods, that are currently under development, use multiple- 
height vector magnetogram data togeth er with the diverge nce-free condition to 
resol ve the azimuthal ambiguity (e.o IWu and Ai Il99fll: Cuperman Li, and 


Semel, 19931 : Li. Cunermari. and Semell. 199.31: Bqulmezaoud and Amari iioo’ol: Li 


Ainari, and Fa n j2007l : Crouch and Barne'Z 20081 : ICrouch. Barnes, and Lekal . l2009l : 

Crouch, 201, ih . Multiple-height vector magnetogram data include information 
about the variation of the magnetic held along the line-of-sight, which can be 
used to approximate the derivatives of the magnetic held in the line-of-sight 
direction. Subsequently, this information can be used to compute the divergence 
of the mag netic held for any position on the solar disk Ifor details see Crouch 
and Barnes, 20081 : ICrouch. Barnes, and Leka 120091 : Crouchl . 2^3l l. However, 
there are several methods available that each make a different assumption about 
how to employ the divergence-free condition to resolve the ambiguity. The dif¬ 
ferences in the implementation details can have a signihcant effect on the results 
produce d: for comparisons and tests using synthetic data see Crouch. Barnes, 
and Leka ( 20091 . henceforth Paper 11 and ICrouc hi (1201.1 henceforth Paper 11). 

All of the methods examined in Paper 1 and Paper 11 use an expression for the 
divergence of the held that involves differentiation of quantities that depend on 
the choice of azimuthal angle. As a result, all heights used to approximate line-of- 
sight derivatives should generally be disambiguated simultaneously. In contrast, 
the purpose of this article is to investigate methods that use an expression for 
the divergence that involves differenti ation of quantities that do no t depend on 
the choice rf,^;Zimuthal^_an£le_[ei^jj_JSkumanic^^n^^eme|j_^9^Semel and 
Skumanich, 19981 : ^kurai and HaginoT 200,3l l. As we will show, this leads to 
an expression for the divergence that can be used to disambiguate each height 
independently. This expression is derived in Section [2] In Section [3] we review 
the synthetic data and metrics that are used to test the perforn iance of the 
variou s disambiguation methods. In Section U] we re-examine the IWu and Ai 
( 1990l l criterion for resolving the azimuthal ambiguity using the expression for 
the divergence derived in Section [2J In Section 0 we test the hybrid method 
that was presented in Paper II, adapted to use the expression for the divergence 
derived in Section [2J In Section [5] we present conclusions. 


ms.tex; 3 November 2009; 3:44; p. 2 






















































































Resolving the Azimuthal Ambiguity in Vector Magnetogram Data 


3 


2. A Modified Expression For the Divergence-Free Condition 


Following ICrouch and Barn^ (2008), Paper I, and Paper II, in this section 
we briefly re-iterate the derivation of the expression for the divergence of the 
magnetic field in terms of observable quantities that is valid for any position 
on the solar disk. We then modify the expression obtained for the divergence 
by replacing derivatives of quantities that are dependent on the choice of az¬ 
imuthal angle with derivatives of quantities that are independent of the choic e 
of azimuthal angle, us i ng an appr oach similar to Skumanich and Serndl (199^, 
Semel and Skumanich ( 1998l l. and Sakurai and Hagiii^ ( 200.81 1 . Subsequently, we 
discuss the consequences for using the modified version of the divergence-free 
condition to resolve the azimuthal ambiguity in vector magnetogram data. 

Throughout this investigation we consider multiple-height, vector-magnetogramj 
data within a field of view that is limited in spatial extent such that a layer of 
constant optical depth can be approximated by the - j/^-plane (where the su¬ 
perscript h denotes heliographic coordinates). At each pixel at each observation 
height, we assume that the image components of the magnetic field vector can 
be measured (denoted with the superscript i); By, and Rl,, where 
is the line-of-sight component and 


i?^ = R^cos^, and By=B^siTi^, (1) 

where B± is the magnitude of the transverse component (perpendicular to the 
line-of-sight) and ^ is the azimuthal angle. It is important to emphasise that, 
when using the linear polarisation of magnetically sensitive spectral lines, the 
azimuthal angle ^ can only be inferred within the range 0 < .^ < 180°, which 
results in the azimuthal ambiguity (i.e., it is not possible to distinguish between 
the two choices of azimuthal angle, ^ and ^ -I- 180°, without additional informa¬ 
tion). The relationship between the heliographic and image components of the 
field is 


B^ = aiiBl, + ai2By + aisB^ , 

By = a2iBl, + a22By + a23Bl, (2) 

B'^ = asiBl, + a32By + assBl, 


where the coefficients a^- are taken from Equation (1) of Gary and HagvardI 


( 1990h and are assumed to be constant within the field of view. 


When the divergence of the magnetic field is expressed in heliographic coor¬ 
dinates, the derivative of the field in the direction perpendicular to the x^~y^- 
plane is required {e.g., dB^/dz^), but this derivative cannot be directly inferred 
from observations except at disk centre. This is because the available methods for 

1 n -PirifVT vnrr +-lTi2i TmT'Tr»+-i r\n r\-F -f /x wi ri rm /x-Fn -R ^1 t /v>h _/ijh T-vloYni^i 


so only in the line-of-sight direction (e.g.. 

Ruiz Cobo and del Toro Iniesta 

1992 

Collados et al. 1994: Metcalf et al. 1995: 

del Toro Iniesta and Ruiz Cobo 

1996 

Liu et al. 1996: Westendoro Plaza et al. 1998 

. 2001: Socas-Navarro, Trujillo 
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Bueno, and Ruiz Cobo 

. 2 OOOI: Elbe et all 20021: Leka and Metcall 20031: Socas- 

Navarro. 2005|. 20071). For convenience, we assume that the variation of the field 


can be inferred as a function of line-of-sight distance z', although we acknowledge 
that in practice the line-of-sight variation of the field may be inferred as a 
function of optical depth r; in such cases, the relationsh ip between r and z' 


can be determined using a solar atmosphere model (e.o ., lMaltby_et_m, [1986 


Vernazza. Avrett. and Loesei , 1981 ; Collados et al. , 1994 Socas-Navarrol . 2007^ . 


Using Equations it can be shown that the relationship between derivatives 
with respect to and is 


df df df 


033 




(3) 


where / is any differentiable function. In practice, all of the derivatives required 
to compute the divergence must be approximated because the field is measured 
at discrete locations: in the horizontal heliographic directions this is because the 
field is sampled discretely over the heliographic plane {e.g., in pixels), whereas, in 
the line-of-sight directi on this is because the magn etic field is generally inferred 
at discrete heights Isee lCrouch and BarnesI ( 2008ll and Paper I for a summary 
of the different methods for inferring the line-of-sight variation of the magnetic 
field). 

Starting with the divergence of the magnetic field expressed in heliographic 
coordinates. Equations © and ([3]) can be used to rewrite the expression in terms 
of observable quantities {i.e., derivatives of the image components of the field 
with respect to and z'), which is valid for any position on the solar disk. 


033 V-B = Da + 033 


where 


dBj 

dz^ ’ 


(4) 


Da = 


dBi dBl dBl^ 

(“11033 - “13031) 

. ^^Bi . . dBl 

((^21(^33 — 0.23(^31) -^^7 + (^^22^33 — 023^^32) -^^7 ■ 


1^^12^133 


(5) 


Each of the terms in Da (Equation dS])) involve differentiation of quantities that 
are sensitive to the choice of azimuthal angle {i.e., B], and By). Eor example, 
in the line-of-sight direction, this means that vector inagnetogram data should 
be disambiguated simultaneously at all heights used to approximate line-of- 
sight derivatives. Paper I and Paper II examined several methods based on 
Equations O and ([5|). The purpose of this article is to investigate methods 
th at are based on a different e x pression. By using an a pproa ch similar to that 
of ISkumanich and Semell (lI996^.ISemel and Skumanichl (Il998ll. and Sakurai and 
Hagino ()2003l l. we modify this expression for the divergence so that differentia¬ 
tion only operates on quantities that are independent of the choice of azimuthal 
angle, such as cos^ sin^^, cos^ sincos (2^), and sin (2^); we prefer to work 
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with these quantities instead of and Bl.By because methods based 

on these perform slightly better in tests (results not shown). For example, we 
can use 


d cos^ ^ 


2cos^ 


i9cos^ 

dx^ 


and 


to show that 


i9sin^^ „ . ,9sin^ 


d cos f sin £ , cos £ sin £ 

= sm£—-^+co&£- 


dx^ 


dx^ 


dx^ 


d cos £ 
dx'^ 


and 


cos.^ 

cos^ 


d cos^ £ 
dx'^ 

9 cos (2,^) 
dx^ 


2 sin^- 
+ sin £ 


d cos £ sin £ 
dx'^ 

9 sin (2^) 
dx'^ ’ 


cos 




d sin^ £ 
dx'^ 


( 6 ) 

(7) 


as^ 

dx^ 


sin^ 

cos^ 


d sin^ £ 
dx^ 

9 sin (2^) 

9^7 


2 cos £ ■ 
— sin^ 


d cos £ sin^ 
9ih 

9 cos(2^) 
dx^ 


. ^ 9 cos^ £ 


( 8 ) 

(9) 


For all of the terms on the right hand sides of Equations © - (0), differentiation 
only operates on quantities that are independent of the choice of azimuthal angle. 
Subsequently, Equations © and © can be used to show that: 


^ dx^ 


= 2Bj_ 


= cos£ 


dcos£ 
dx^ ^ 
9Rj. 
dx^ 


2 cos £ 


dBi 


B^ 


dx^ 

9 cos(2^) 


dx^ 


+ B± sin £ 


9 sin {2£) 
dx^ 


( 10 ) 


and 


dBl 

0 _ y 

dx^ 


= 2B^ 


= sin^ 


9sin^ 
dx^ 
dB± 
dx^ 


2 sin £ 


dB^ 


-B^ 


dx'^ 

9 cos(2^) 
9a;*' 


+ B± cos £ 


9sin(2g) 

9a;*' 


( 11 ) 


Again, for all of the terms on the right hand sides of Equations (fTn|) and (ITT]) , 
differentiation only operates on quantities that are independent of the choice of 
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azimuthal angle. Analogous expressions can be derived for the derivatives of B^. 
and By with respect to both and z‘. Using Equations m and m and their 
counterparts, we can re-write the divergence-free condition as 


where 


2033 V-B 


hs sin ^ + he cos ^ -I- 2033 



( 12 ) 


and 


^ „ 5 sin (20 

hs — 031-DU- TT-: -1-032 


dz^ 

+ (011033 — 013031)5^ 

-I- (012033 — 013032) 


cos (20 


dz 
9 sin (20 
dx^ 


dz' 


+ (021O33 — 023O31) B ± 


p dBjL _ 9cos(20 

dx^ ^ dx^ 

9 sin (20 


dy^ 


(022033 — 023 O 32 ) 


2 ^^ a cos (20 


dy^ 


dy^ 


(13) 


he — 


O 31 


dB± 


B^ 


dz' 

(011O33 — 013031) 


a cos (20 

a^i 

dB± 


0-32B± 


a sin (20 


dz' 


+ (012O33 — 013032) B± 


(021033 — 023O31) 


-I- ( 022 O 33 — 023 O 32 ) B± 


dx^ ' 
a sin (20 

a sin (20 


a cos (20 
dx^ 


a cos (20 


dy^ 


dy^ 


(14) 


The quantities hg and he are independent of the choice of azimuthal angle. 
Moreover, in hg and he, differentiation only operates on quantities that are inde¬ 
pendent of the choice of azimuthal angle. This means that the approximation for 
the divergence of the magnetic field at one location does not need to depend on 
the choice of azimuthal angles in neighbouring locations (if the various quantities 
are evaluated at the centre of each pixel at each height). Therefore, it is possible 
to develop disambiguation methods based on Equations (fT^ - (ITTl) such that 
the choice of azimuthal angle at one location does not depend on the choice of 
azimuthal angle at other locations (in any direction). This can have a significant 
impact on methods for resolving the azimuthal ambiguity with the divergence- 
free condition. Eor example, in the line-of-sight direction, using a method based 
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on Equations (HU) - (m, each height can be disambiguated independently; this 
is different to the methods examined in Paper I and Paper II, which were based 
on Equations m and ([5|). It should be noted that line-of-sight derivatives of 
cos (2^), sin (2^), B±, and are generally required to compute the divergence 
using Equations (na) - (HI, but these do not depend on the choice of azimuthal 
angle. 


3. Synthetic Data and Performance Metrics 


To test the performance of the various algorithms based on the divergence-free 
condition, as described by Equations (Ha - (Ha, we use the same synthetic data 
that were used in Paper I and Paper II for which the correct configuration of 
azimuthal angles is known and measurements of the magnetic field are available 
at two heights. These synthetic d atasets are based on the s ame magnet i c field 
configurations that were used bv iMetcalf et all ( 2006l l and iLeka et o,l. i 20091 1 
and, therefore, allow direct comparisons with the single-height methods tested 
in those articles. Three classes of synthetic data are considered: 

i) A snapshot from the simulations of iFan and GibsonI ( 2004 ) of a twisted 
magnetic flux tube emerging into an overlying potential field arcade. This 
magnetogram is situated at disk cen tre and does not inc lude the effects of 
noise. Further details can be found in iMetcalf et al. (200^ and Paper I. This 
synthetic dataset was provided by Yuhong Fan. 

ii) A multipole, linear-force free magnetic field constructed from a set of point 
sources located on a plane below the solar surface, located away from disk 
centre, centred at (9°S, 36°E). Three cases are constructed with different levels 
of noise added to the polarisation spectra: one no-noise case and two noise- 
added cases. Further details are provided in Leka et al. 1 2009h and Paper II. 
These synthetic datasets were provided by K.D. Leka and Graham Barnes. 

iii) A potential magnetic field with fine-scale structure with properties similar to 
penumbrae and plage, located at disk centre, and sampled on a grid with a 
pixel size of 0.03". To simulate the effects of imperfect instrumental spatial 
resolution three magnetograms are constructed with the polarisation spectra 
spatially binned by factors of 5, 10, and 30, resulting in effective pixel sizes 
of 0.15", 0.3", and 0.9", respectively. Some of the fine-scale structures that 
are fully resolved at 0.03" are not fully resolved at lower spati a l reso lution. 
For some additional discussion of this test case see Leka et al\ (120091') . Leka 
and Barnes ( 2012l L GeorgouliH 2012l L Leka et all ()2012l l and Paper II. These 
synthetic datasets were also provided by K.D. Leka and Graham Barnes. 

We emphasise that obtaining measurements of the magnetic held at two 
heights for solar observational data is a challenge that is beyond the scope of this 
investigation. In principle, the li ne-of-sight variation of the magnetic he ld can be 
inferredJromobservationsfSiSiJfiS^^S^^^S^^sI^i^sissSii^^ LColl ados 


et al.. 199417 Metcalf et al. 1995 Hdel Toro Iniesta and Ruiz Gobol. 19961: Liu 


et al, 19961 : IWestendorp Plaza et all igpsT 200ll: Socas-Nayarro Trujillo Bueno 
and Ruiz Cobo, 2000l : Eibe et aTT 20021 : Leka and Metcalf . 20031 : Socas-Navarrol . 
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l2005l . [200^ . However, further research is required to develop methods for reliably 
obtaining maps of the solar magnetic field vector at multiple heights; this is 
especially true if one of the observation heights is located in the chromosphere 
(where interpretation difficulties are involved). 

There are other sources of uncertainty that are not modelled by the synthetic 
data used in this investigation that may be significant for solar observational 
data: i) it is assumed that the location of the observation height is known exactly 
(be., the effects of limited spatial resolution in the line-of-sight direction are 
ignored); ii) the magnetic field is assumed to be measured at constant geometrical 
height (constant optical depth may be more appropriate); and iii) the physical 
distance between the two observation heights is assumed to be constant (but this 
may vary over the field of view). It is important to test the influence of these 
sources of uncertainty before the methods presented in this article are applied 
to solar observational data. 

To measure the performance of the algorithms we use: Alarea, the fraction of 
pixels correctly disambiguated, and A4 the fraction of transverse magnetic 
held correctly disambiguated above a threshold T] we present results for T = 
100 G (gauss) and T = 500 G. 


4. The Wu and Ai (1990) Criterion 


The idea behind the IWu and Ai 1 1990l l criterion for resolving the azimuthal 
ambiguity is to express the equation for t h e div e rgence-free condition as an 
inequality (e.q., Cunerman. Li, and SemdI. 1^93 Li. Cuperman. and Sem5l . 


Il993l: III Amari. and Fanl . 120071 : iGrouch and Barnesl . l2008l : Paper I: Paper III . 

For example, assuming V-B = 0, Equation (IT^ can be used to show that 


033 -^ {hsSin^ + hcCOS^} = -2 1 “ 33 -^ 


< 0 . 


(15) 


Thus, given the sign of dB^^jdz\ Equation (ITSl) may be used to resolve the 
ambiguity, by choosing the realisation of the azimuthal angle that produces the 
correct sign for the quantity hg sin^ + he cos^. The magnitude of dB'^/dz' is not 
required; however, regarding various line-of-sight derivatives, it should be noted 
that the sign and magnitude of the line-of-sight derivatives of B±, cos (2^), and 
sin (2^) are required to implement Equation (1151) . except for locations where 
031 = 0 and 032 = 0 . 

We implement the Wu and 'Ail dlbblt criterion using Equation m as follows. 
Hereafter, pixel labels in the x‘-, y‘-, and z‘-directions are given by i, j, and k 
respectively, and the number of pixels in the x*- and j/'-directions are given by Ux 
and Uy, respectively. For pixels with 1 < i < Ux and 1 < j < Uy, the horizontal 
heliographic derivatives d/dx^ and d/dy^ at pixel (i, j, k) are approximated with 
three-point finite differences using measurements of the magnetic field at pixels: 
{i,j, k), (f -I-1, _), k), and {i,j + l,k). At the boundaries where i = rix and j = riy, 
the approximations for the horizontal heliographic derivatives are modified to 
use only pixels within the field of view. At pixel (f, j, k) line-of-sight derivatives 
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djdz^ are approximated with forward finite differences using measurements of the 
magnetic field from two heights, at pixels {i,j, k) and (i,j, k + 1). Subsequently, 
at each pixel the various terms in Equation (1151) are computed. If the inequality 
is satisfied the choice of azimuthal angle is not changed. Otherwise, the choice 
of azimuthal angle is switched. Three potentially advantageous features of this 
approach are i) the disambiguation result at one location does not depend on 
the choice of azimuthal angles in neighbouring locations; ii) it is non-iterative, 
meaning that each pixel only needs to be visited once during disambiguation; 
and iii) it is non-sequential, in that the results do not depend on the order in 
which pixels are visited. 

Some limitations of the Wu and Ai 1 1990f ) criterion implemented with Equa¬ 
tions © and were demonstrated in Paper I. These include i) the disambigua¬ 
tion results depend on the smoothness of the initial configuration of azimuthal 
angles; ii) the disambiguation results depend on the order in which pixels are 
visited; and iii) the approximation for Da in Equation ([S]) can have the same 
sign for both the correct and the incorrect choices of azimuthal angle (when 
deriva tives are approxim ated from discrete measurements). The implementation 
of the IWu and Ail (|l990ri criterion based on Equation (na, as described above, 
does not suffer from these problems. This is because the derivatives in and he 
(see Equations (USD and (iMll . respectively) only depend on quantities that are 
independent of the choice of azimuthal angle. Furthermore, because hs and he 
are independent of the choice of azimuthal angle, the quantity hs sin ^ + he cos ^ 
changes only in sign (not magnitude) when the choice of azimuthal angle is 
switched at a given pixel at a given height. However, the divergence of the mag¬ 
netic field is generally not exactly zero when derivatives are approximated from 
discrete measurements and, therefore, some erroneous solutions are expected. 



Figure 1. Results for the various methods applied to the twisted flux tube and arcade. Areas 
with the correct/incorrect azimuthal angle are black/white. Positive/negative vertical magnetic 
flux is indicated by red/blue contours at 100, 1000, 2000 gauss (G). The magnetic neutral line 
is indicated by the green contour, (a) Results from the Wu and Ai (1990) criterion, (b) Results 
from the global minimisation method, which minimises Fg (see Equation ([Teli'i. (c) Results 
from the hybrid method, which first minimises Eg as in (b) (see Equation l llOIl h then applies 
a smoothing algorithm to pixels below a threshold transverse field strength (as described in 
Section 0. 


For the flux tube and arcade (Figure [T](a) and Table [T|), the IWu and~Ai 
( 199nl l criterion implemented with Equation (11511 performs significantly better 
than the implementation used in Paper I (with Equations © and ©) in that 
the assumption made is correct over f or most of the fiel d of view. It should be 
noted that the assumption made by the Wu and Ail ( 19901 1 criterion implemented 
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with Equation (1151) is correct/incorrect when the correct/incorrect realisation 
of the azimuthal angle is selected; that is, the sign of the quantity /igSin^ + 
he cos 5 is correct for the correct choice of azimuthal angl e. Likewise, f or the 
multipole field without noise (Figure [5] and Tabled, the IWu andAil (|l99fll 'l 
criterion implemented with Equation (HH) performs significantly better than the 
implementation that uses Equations and as demonstrated in Paper I and 
Paper n^loweverj_as^heJevelof^ioisejncreaseSjJhejesults^roducedb^th^Wu 
and Ai (|l99(i l criterion implemented with Equation (fTKI) degrade dramatically 
(Figure [2] and Table [2]). Undesirable results are also produced for the cases with 
limited spatial resolution (Figure [3] and Table [3]). 


Table 1. Performance metrics for the various methods applied to the twisted 
flux tube and arcade. 



■Ad area 

AIbj^>ioo g 

MB _ e >500 G 

Wu and Ai fl990') 

0.98 

0.98 

0.97 

Global minimisation, |V-B| S 
(Equation Qbjl) 

1.00 

1.00 

1.00 

Hybrid 

1.00 

1.00 

1.00 


We have investigated an approach that assumes that the divergence of the 
magnetic field is minimised at each pixel, approximating derivatives as described 
above. This approach is non-iterative, non-sequential, and the solution at one 
location does not depend on the choice of azimuthal angles in neighbouring 
locations. However, when implemented with Equations (1X31) ~ (I14L it c an be 
shown that such an approach produces identical results to the lWu and Ail ( 199Cll l 
criterion implemented with Equation (1151) . 


5. The Hybrid Method 

In the previous section, we found that methods for which the solution at one 
location does not depend on the choice of azimuthal angles in neighbouring 
locations tend to produce undesirable results when either noise or unresolved 
structure are present in the magnetogram data. In this section we adapt the two- 
step, hybrid method that was presented in Paper II, for which the solution at one 
location does depend on the choice of azimuthal angles in neighbouring locations. 
The first step of this approach involves the global minimisation of a combination 
of the approximation for the divergence and a smoothness constraint. In the 
second step, pixels with a weaker transverse component of the magnetic field are 
revisited with a smoothing algorithm. 

5.1. Step One: The Global Minimisation Method 

For the first step of the hybrid method, we assume that the correct configuration 
of azimuthal angles corresponds to the minimum of 
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Figure 2. Results for the various methods applied to the multipole field with different 
degrees of noise. Areas with the correct/incorrect ambiguity resolution are black/white. Posi¬ 
tive/negative vertical magnetic flux for the noise-free case is indicated by red/blue contours at 
50, 700, 1500 and 2800 G. The magnetic neutral line is indicated by the green contour. (Left to 
right) No-noise, low-noise and high-noise cases. (Top row) Results from the Wu and Ai (1990) 
criterion. (Middle row) Results from the global minimisation method, which minimises Es (see 
Equation (I16II ). (Bottom row) Results from the hybrid method, which first minimises Es as 
in the middle row (see Equation ifTOf). then applies a smoothing algorithm to pixels below a 
threshold transverse field strength (as described in Section!^. 


Ux riy 

= + , (16) 

i=l 3 = 1 

where (V-B)ij_fc is the approximation for the divergence of the magnetic field 
at pixel implemented using Equations (I12|) - (THll with derivatives ap¬ 

proximated as described in Section|4j As is a dimensionless parameter, and Sij^k 
is given by 
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Table 2. Performance metrics for the various methods applied to the multipole field with different degrees! 
of noise. I 


Metric: 

Noise level: 

None 

Ad area 

Low 

High 

M 

None 

S ^>100 

Low 

G 

High 

M 

None 

Bj ^>500 

Low 

G 

High 

Wu and Ai fl990^ 

0.92 

0.61 

0.54 

0.95 

0.71 

0.57 

0.97 

0.80 

0.61 

Global minimisation, | V-B| + S 
(Equation lllfilll 

1.00 

0.99 

0.94 

1.00 

1.00 

0.98 

1.00 

1.00 

1.00 

Hybrid 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


Si j k — 


121 1/2 


It' • 2 • 21"^ 

+ +[Bl{i,j,k) - B]y{i + l,j,k)\ I 

1 r • 9 • 91 1 /2 


where di is the distance between the centres of the pixels at {i,j, k) and {i + 1, j, fc),| 
and d 2 is the distance between the centres of the pixels at {i,j, k) and (i, j + 1, fc). 

At the boundaries of the field of view, where i = Ux and j = Uy, the definition 
of Si^j^k is modified to use only pixels within the field of view. The purpose of 
Si^j^k in Equation m is to minimise the difference between the magnetic field in 
neighbouring pixels. Thus, at a given pixel the first term in Equation (uni) does 
not depend on choice of azimuthal angles in neighbouring pixels whereas the 
second term does (in the horizontal heliographic directions but not in the line- 
of-sight direction). Because the definition of Sij^k (see EauationfTTl) involves only 
measurements of the magnetic field at a single height k, each height can be dis¬ 
ambiguated independently. Consequently, this approach is less computationally 
expensive than the corresponding approach presented in Paper II, which disam¬ 
biguates all heights simultaneously because it was based on Equations and 
®. For example, to disambiguate two heights using this approach takes about 
half the compute time of the corresponding approach presented in Paper II. 

To find the configuration of azimuthal a ngles that corresponds to the mini¬ 
mum of Eg we use simulated annealiu! 


mum ot fig we use smuUatedanneaimgj_ej^ 
Gelatt, and Vecchi, 19831 : Press et al\ . 1992l l. 


Metropolis et al\. 19531: Kirkpatrick. 


The simulated annealing algorithm 
is implemented exactly as described in Paper II. We disambiguate each magne¬ 
togram twenty times with different sequences of random numbers and select the 
solution that yields the lowest value of E^. The value of the parameter Ag does 
affect the solution obtained; for simplicity, we set Ag = 1 for all cases. Results 
from this first step of the hybrid method are presented in Figure[TIb), Figure[2] 
(middle row), and Figure [3] (middle row), and Tables (T] [31 and [3] 

For the flux tube and arcade (Figure [T](b) and Table |T|) many of the solutions 
retrieved by the global minimisation method have a value of Eg less than that 
for the correct solution. This indicates that the assumption that the correct 
configuration of azimuthal angles corresponds to the global minimum of Eg 
(Equation p6ll l is violated in these cases. However, for the solution with the 
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Figure 3. Results for the various methods applied to the limited resolution cases. Areas with 
the correct/incorrect ambiguity resolution are black/white. The green contour is the magnetic 
neutral line. (Left to right) 0.15'', 0.3'' and 0.9" cases. (Top row) Results from the Wu and Ai 
(1990) criterion. (Middle row) Results from the global minimisation method, which minimises 
Es (see Equation {Hi). (Bottom row) Results from the hybrid method, which first minimises 
Es as in the middle row (see Equation 1 II 6 I 1 I. then applies a smoothing algorithm to pixels 
below a threshold transverse field strength (as described in Section [5ll . 


lowest value of Eg, the incorrect azimuthal angle is only obtained in a small 
number of pixels. 

For the multipole field without noise (Figure [5] and Table [2]) the global 
minimisation method retrieves the correct solution at every pixel. However, as 
the noise level increases, the performance of the global minimisation method 
decreases. This is most evident in regions where the magnitude for the transverse 
component of the field is relatively weak and the signal tends to be dominated 
by noise {e.g., Ki 0; x' Ki 300; and at x' ss 170 and y' Ri 50). For both noise- 
added cases, the solutions retrieved by the global minimisation method have a 
value of Eg less than that for the correct solution. Nevertheless, the quality of 
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Table 3. Performance metrics for the various methods applied to the limited resolution cases. | 


Metric: 

Pixel size: 

0.15" 

Ad area 

0.3'' 

0.9" 

Mb 

0.15" 

_L>100 

0.3" 

G 

0.9" 

Mb 

0.15" 

_L>500 

0.3" 

G 

0.9" 

Wu and Ai I'lQQO) 

0.96 

0.93 

0.84 

0.96 

0.93 

0.82 

0.96 

0.93 

0.81 

Global minimisation, | V-B| + S 
(Equation 

0.99 

0.99 

0.99 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 

Hybrid 

0.99 

0.99 

0.99 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 


the solutions retrieved by the global minimisation method that minimises Eg is 
very similar to that obtained by the corresponding method from Paper II. 

For all of the cases with limited spatial resolution (Figure |3] and Table [3|) 
we again find that the assumption made by the global minimisation method is 
violated (be., the solutions retrieved have values of Eg less than that for the 
correct solution). However, the quality of the solutions retrieved by the global 
minimisation method is generally comparable to the corresponding method from 
Paper II. The approach implemented here tends to do better than the corre¬ 
sponding method from Paper II in the plage region (at the centre, right of the 
field over field), and slightly worse around the flux concentration at the top, 
right of the field of view and between the two flux concentrations in the bottom, 
centre of the field of view. All of these are regions with significant unresolved 
structure (see Fi gure 2 of Paper II) and pose a challenge for all methods (see 
Leka et al. . 2009l l. 


5.2. Step Two: A Smoothing Algorithm 


Evi dently, the g l obal m inimisation method produces much better results than 
the Wu and Ai ( IQQOh criterion. However, it produces undesirable results in 
regions where the transverse component of the magnetic field is strongly affected 
by photon noise. Consequently, in the second step of the hybrid method we 
apply a smoothing algorithm that revisits pixels in regions where the transverse 
component of the field is weak. The smoothing algorithm does not use the 
divergence-free condition. In contrast, it aims to minimise the difference between 
the magnetic field in a small neighbourhood of pixels. The smoothing algorithm 
is identical to the method described in Section 5.4 of Paper H, except that the 
definition of the neighbourhood is modified to include only pixels from a single 
height. This allows each height to be disambiguated independently. For the no¬ 
noise, low-noise, and high-noise cases we revisit pixels with a transverse field 
strength of 100 G, 150 G, and 400 G, as in Paper H. 

For the flux tube and arcade (Figure [IJc) and Tabled]), the smoothing algo¬ 
rithm corrects the solutions in the small number of pixels that were incorrectly 
resolved by the global minimisation method. The resulting solution has the cor¬ 
rect azimuthal angle for all pixels. For the multipole field without noise (Figure^ 
and Tabled]), the smoothing algorithm does not change the solution in any pixels. 
For the noise added cases the smoothing algorithm substantially improves the 
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solutions in the regions that are revisited. For the cases with limited spatial 
resolution (FigureOand Table[3]), the smoothing algorithm changes the solution 
in a small number of pixels in each case, but does not substantially change the 
overall quality of the solutions. Generally speaking, the final solutions produced 
by the two-step, hybrid method implemented here are very similar in quality to 
those produced by both the corres ponding approach from Paper II and the better 
performing methods examined in Leka et al. ( 2009ll . The smoothing algorithm 
works well for the synthetic datasets used in this article, which are quite limited 
in spatial extent; however, we note that it may not be the optimal method for 
use over extended areas far from regions where the solution is obtained from the 
global minimisation method. 


6. Conclusions 


We continue the investigation of how to use the divergence-free condition to 
resolve the azimuthal ambiguity in vector magnetogram data. In previous work 
(Paper I and Paper II), all of the methods examined used an expression for the 
divergence of the magnetic field that involves differentiation of quantities that 
depend on the choice of azimuthal angle. As a result, all heights used to approxi¬ 
mate line-of-sight derivatives should generally be disambiguated simultaneously. 
In contrast, in this article we investigate a set of methods that use an expression 
for the divergence that involves differe ntiation of quantities that do not depend 
on the choi^_o£azimutha^^ngle_^e;^;jlSkumanicl^^n^_Sem^ 199^Semel and 
Skumanich, 19981 : ^kurai and Haginol . boOSll . This results in an expression for 
the divergence that can be used to disambiguate each height independently. 

We test two methods that each make different assumptions a bout how to 
use th e divergence-free condition to resolve the ambiguity: i) The IWu and Ai 
( 199nh criterion, whi ch expresses the equation for the divergence-free condition 
as an inequality {e.q., Cuperman, Li, and Semel, 1993: lLi. Cunerman. and SemdI . 


1993 : Li. Amari. and Fanl . 20071 : Crouch and Barnes . 2008 : Paper I; Paper II); 
and ii) The two-step, hybrid method from Paper II, adapted to disambiguate 
each height independently, which substantially improves its computational ef¬ 
ficiency. The first step of the hybrid method involves the minimisation of a 
combination of the approximation for the divergence and a smoothness con¬ 
straint. The second step revisits pixels with a weaker transverse component of 
the magnetic field with a smoothing algorithm. To test these methods we use 
three different types of synthetic data: i) error-free synthetic data; ii) synthetic 
data with noise to simulate Poisson photon noise in the observed polarization 
spectra; and iii) synthetic data that includes a spatial binning to simulate the 
effects of limited instrumental spatial resolution in t he directions perp endicular 
to the line-of-sight. For error-free synthetic data the Wu and aH (1990) criterion 
implemented in this article performs better than the implementation in Paper I 
a nd Paper II. However, we again find that the results produced bv the Wu and 


AiH 


criterion are very sensitive to both noise and un resolved structure . 


We also find again that the hybrid method outperforms the lWu and Ai ( 1990l ') 
criterion. Furthermore, the quality of the results retrieved by the hybrid method 
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implemented in this article are generally very similar to those retrieved by the 
corresponding method in Paper II, yet the compute time taken to produce the 
results is substantially reduced. The software for the global minimisation method 
and the hybrid method is available at http://www.cora.nwra.com/~ash/ambigl. 
tar.gz. 
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